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ABSTRACT 

We present the star formation rate (SFR) and starburst fraction (SBF) for a sample of field galaxies from 
the ICES intermediate-redshift cluster survey. We use [O II] and Spitzer 24 fim fluxes to measure SFRs, and 
24 fj,m fluxes and HS absorption to measure of SBFs, for both our sample and a present-epoch field sample 
from the Sloan Digital Sky Survey (SDSS) and Spitzer Wide-area Infrared Extragalactic (SWIRE) survey. We 
find a precipitous decline in the SFR since z = L in agreement with other studies, as well as a corresponding 
rapid decline in the fraction of galaxies undergoing long-duration moderate-amplitude starbursts. We suggest 
that the change in both the rate and mode of star formation could result from the strong decrease since z = 1 of 
gas available for star formation. 

Subject headings: galaxies: evolution — galaxies: starburst — galaxies: stellar content 



1. INTRODUCTION 

In recent years, much attention has been focused on the 
early history of star formation, in particular, the rise of the 
SFR to its peak around z - 2 (e.g., Fafl, Chai'lot, & Pei 1996, 
Madau et al. 1996, Giavalisco et al. 2004, Perez-Gonzalez 
et al. 2008, Bouwens et al. 2008). Such a rapid rise is not sur- 
prising considering the dynamical timescales of galaxy-sized 
structures and the lifetimes of stars. More surprising, perhaps, 
is the precipitous decline in the global SFR since z ~ 1 — 
more than a factor of 10 lower than its peak value, and falling 
fast (e.g., Lilly et al. 1996, Schiminovich et al. 2005, Hop- 
kins & Beacom2006, Villar et al. 2008). This is remarkable 
not only for its rapidity but also because it appears to mark our 
epoch as the beginning of the end of galactic star formation. 
It is now clear that the unexpected prevalence of starform- 
ing galaxies in rich galaxy clusters discovered by Butcher & 
Oemler (1978), long regarded as a cluster phenomenon, is 
universal, as starforming galaxies — great and small — sur- 
render youthful vigor and fade towards oblivion in only a few 
billion years. 

Using the wide field of the Inamori-Magellan Areal Camera 
and Spectrograph (IMACS) on Magellan-Baade, the MACS 
Cluster Building Survey (ICBS) is focused on the study of 
galaxy infall and evolution from R 5 Mpc into cluster 
cores. Because the projected density of cluster/supercluster 
members is low at such large radii, our near-complete sam- 
ples necessarily include ^1000 "field" galaxies at redshift 
0.2 < z < 0.8 per survey field. This gives us an opportunity to 
compare galaxy evolution in clusters with the field over this 
epoch. From these data we report in this Letter on the signifi- 
cant decline in SFR and SBF since z ^ 1 . 

2. DATA 

The data discussed here come from 4 fields that contain rich 
galaxy clusters at z = 0.33,0.38,0.42, &0.55. The IMACS f/2 
spectra have an observed-frame resolution of 10 A full-width- 
half-max with a typical S/N ~ 20-30 in the continuum per 
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resolution element. Spectral coverage varies, but almost all 
cover [O II] and H/3 emission, and HS absorption, our opti- 
cal diagnostics of star formation; a fraction cover Ha as well. 
In each 28'-diameter IMACS field we have observed 65% of 
the galaxies that are brighter than r ^ 22.5, obtaining ade- 
quate spectra of 8 1 % of these. Measurement of spectral fea- 
tures followed procedures described in Dressier et al. (2004, 
hereinafter D04). For two fields we have confusion-limited 
Spitzer Multiband Imaging Photometer (MIPS) 24 fim im- 
ages (Guest Observer program 40387) that cover nearly the 
entire IMACS fields. Details of the data, data reduction, 
and analysis of the field sample are described in Oemler 
et al. (2009b, hereinafter Oem09b). 

In the following analysis we include all galaxies between 
0.10 < z < 0.70 with absolute AB magnitudes at 4400 A 
brighter than M44 + 1 .0, except for those with redshifts within 
±3(7 of each targeted cluster M44 has been determined to 
evolve with redshift as M^^ = -20.00 -z (Oem09b); the Hmit 
of + 1 .0 is about equivalent to r = 22.5 at z = 0.70. There 
are 1 144 objects in this sample. Galaxies are given weights 
proportional to the inverse of the incompleteness of the data at 
each apparent magnitude. To provide a low-redshift point we 
use two samples of SDSS galaxies with redshifts 0.04 < z < 
0.08. To compare to infrared-derived properties of the ICBS 
sample, we use a sample of 385 SDSS galaxies within the 
SWIRE fields (Lonsdale et al. 2003). For comparison with 
the optically-derived properties of our sample, we use this 
SDSS/SWIRE sample and add 690 SDSS galaxies near the 
North Galactic Pole. We assume a concordance cosmology 
with ftmaiter = 0.27 and Ho = 71km s"' Mpc"'. 

3. MEASUREMENTS OF SFR SINCE Z = 0.7 

For galaxies with Spitzer-MIPS coverage we calculate 
SFRs in a way similar to Perez-Gonzalez et al. (2006), us- 
ing the 24 /im flux to estimate the absorbed UV flux and 
[O II] to estimate the escaped UV flux, but including the k- 
corrections of Rieke et al. (2008) and a SFR scale based on 
a Salpeter initial-mass-function (IMF, vid. Oem09b). Donley 
et al. (2008) report that 10-15% of 24 /im sources at these 
flux densities are AGN-dominated. Given the strong redshift 
evolution of the AGN luminosity function, the contamina- 
tion in our z < 0.7 sample should be should be considerably 
lower. Figure 1 (bottom panel) shows the median SFR per 



unit L/L* , which we shall refer to as the specific SFR, as 
a function of redshift. We choose the median specific SFR 
over the more commonly reported mean, because it is pro- 
vides a more stable and typical value — particularly at higher 
redshifts, where the mean is dominated by a small number 
of very luminous objects. (In the 0.60 < z < 0.70 interval, 
only 1.6% of the objects contribute 56% of the total star for- 
mation.) We also show, for comparison, values of the me- 
dian SFR/(L/L*), derived from the mid-IR data of Damen 
et al. (2009), and from the Perez-Gonzalez et al. (2008) study 
of the buildup of stellar mass (from near-IR measurements), 
as reinterpreted by Damen et al. and converted by us to a 
Salpeter IMF scale. To convert the Damen et al. values of 
SFR per unit mass to SFR/ (L/L*) we take a mean value of 
M/L at 4400 A of 3, which is the ratio of the cosmic mass 
density (Perez-Gonzalez, et al. 2008) to the cosmic luminos- 
ity density (Blanton et al. 2003). Although cosmic means, 
both are dominated by galaxies in the luminosity and mass 
ranges of our sample, and thus appropriate for our use. With 
M/L = 3, a sample with our absolute magnitude limit has a 
mean mass of 5a;1O"^M0. 
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Fig. 1 . (Bottom) Median Specific SFR, in solar masses per year per L* luminos- 
ity. Black filled circles show values for our sample derived from Spitzer-MIPS 24 /^m 
fluxes. Red stars are data from Damen et al. (2009), and red open circles are data from 
Perez-Gonzalez et al. (2008) as reinterpreted by Damen et al.(2009). The median SFR 
appears to be declining more rapidly towards the present epoch. (Top) Median values 
(solid circles) and 90th percentile values (open circles) of [O UJfor the z > 0.2 ICBS 
field sample, binned to achieve comparable of numbers of [O II] -detected galaxies. The 
single point at ^ = 0.06 comes SDSS galaxies with optical spectroscopy (see text). 

At z > 0.4 the median values of SFR increase more slowly 
than the mean SFR, as reported by Damen et al. (2009), Mar- 
tin et al. (2007), Zheng et al. (2007), but are otherwise quali- 
tatively similar However, what has not been apparent from 
previous observations is that the fall of the median SFR ap- 
pears to be accelerating in recent epochs. 

Optical emission lines [O II] , H(3 , and Ha are less reli- 
able SFR indicators than 24 ^m flux because of the very large 
and variable extinction which blankets the HII regions of star- 
forming galaxies. However, they do provide at least a qualita- 
tive measure of current star formation and, unlike our Spitzer- 
MIPS data, are available for all four ICBS fields. Therefore, 



we plot in the top panel of Figure 1 the 50th and 90th per- 
centiles of the distribution of equivalent width EW([0 II]) as 
a function of redshift. The shapes of these distributions are 
qualitatively similar to the infrared-derived SFR, however, 
[O II] strength increases more slowly with redshift because 
of the well-known effect that galaxies with higher SFRs also 
have higher dust extinction. 

4. WERE STARBURSTS THE NORMAL MODE OF STAR 
FORMATION BEFORE Z = 1? 

It is now well established that starburst and post-starburst 
galaxies are abundant in intermediate-redshift clusters (Pog- 
gianti et al. 1999 (hereinafter P99), Oemler et al. 2009a), how- 
ever, less is known about the prevalence of starbursts among 
field galaxies at these epochs. Starburst indicators include ex- 
ceptionally strong Balmer absorption, intense optical emis- 
sion, and excess 24 jim flux. Strong Balmer absorption in 
the integrated light of galaxies signals a rapid decline in the 
SFR, because light from A stars persists after the blue contin- 
uum from O and B stars, and the emission from HII regions 
— both of which dilute the Balmer absorption lines — be- 
gin to fade or disappear altogether (P99). Although strong 
Balmer absorption accompanies even the simple truncation of 
star formation, as O and B stars and their HII regions evolve 
away, stellar EW(}\5 ) > 5A occurs only in the aftermath of a 
significant rise in the SFR above its past average, i.e., a burst 
(Dressier & Gunn 1983, Couch & Sharpies 1987, P99). Us- 
ing this criterion, P99, D04 and Oemler et al. (2009a) found 
that the SBF in the field at z ~ 0.4 is lower than in clusters, 
but significantly higher than it is in the field today; however, 
the data sets used were quite small. More recently, Poggianti 
et al. (2009) find a large incidence of dusty starburst galax- 
ies for 0.4 < z < 0.8 in all environments, but particularly 
in groups. Using an entirely different approach. Bell et al. 
(2005) find a substantial fraction of massive field galaxies at 
z ^ 0.7 with SFRs much higher than their long-term aver- 
ages. On the other hand, Noeske et al. (2007) argue that the 
narrow width of the SFR/mass versus mass relation among 
field galaxies in the Groth strip of DEEP2 data is inconsistent 
with a large fraction of strong starbursts. 

Determining the starburst frequency — SBF — as a func- 
tion of cosmic epoch and environment is important both for 
understanding the nature and cause of starbursts, and their role 
in galaxy evolution. For example, starbursts provide a quick 
means of exhausting the gas supply in a galaxy. The dense en- 
vironment of rich clusters provides a variety of mechanisms 
for producing starbursts, including ram-pressure from the in- 
tergalactic medium, tidal encounters between galaxies, and 
mergers and accretion. Some of these are also viable in the 
group environment, to more or less effect, but ram pressure, 
for example, is likely to be unimportant. Some present-epoch, 
truly isolated galaxies appear to have experienced starbursts 
as well, which might point to some instability in the disks of 
star forming galaxies, or accretion of a small satellite, as other 
possible causes. 

Figure 2 compares results from two different methods of 
determining the SBF. The top panel shows the fraction of all 
galaxies in our field sample that have EW(H(5 ) of 2.5A or 
greater than expected for a normal, non-bursting galaxy, fol- 
lowing the method of D04, but using an improved EW(H(5) 
versus EW([OII]) relation (q.v. Oem09). Following Bell 
et al. (2005), we also calculate the SBF from the Spitzer- 
derived SFRs, by defining a starburst galaxy as one whose 
observed SFR is significantly higher than its long-term aver- 



age rate, which should be, ignoring mass loss during stellar 
evolution. 



SFR 



past average 



= Mgal/tgal = (Lgal X M/L)/tgai 



(1) 



where t^„i is the length of time that a galaxy has been forming 
stars, and Mgai, Lgai, and M/L are, the galaxy's stellar mass, 
luminosity, and stellar mass-to-light ratio. We define tgai to 
be time elasped from redshift z = 4 (a reasonable assumption 
for the start of star formation) to the observed redshift, and 
assume that Mg^i /Lgai = 3, where Lgai is calculated at 4400 A. 
We identify a galaxy as a starburst if SFRohs/ SFR past-average > 
3. We plot in Figure 2 the fraction of galaxies which meet this 
criterion and perform the same analysis for the observations 
of Noeske et al. (2007) and Bell et al. (2005). For these we 
include objects with Mgai > 1.6 x lO'^M©, which for M/L = 
3, is equivalent to our luminosity Umit of + 1 .0. 



same method as we and Bell et al. (2005) have used, is in good 
agreement with our own. 

Finally, we show another way of assessing the SBF(z), us- 
ing the composite spectra approach developed in D04. Here 
we add, weighted by luminosity, all ICBS spectra in each of 5 
redshift slices 0.20 < z < 0.70, and for a present-epoch sam- 
ple of SWIRE galaxies. We plot in Figure 3 EW([0 II] , US) 
measured from these composite spectra, along with expected 
values for a population of (mostly) continuously-starforming 
galaxies at the present epoch. Figure 3 includes expected 
values calculated by D04 from the 2dF & CN0C2 field 
surveys (see D04) and also new determinations using the 
SWIRE data, an improvement because we now have individ- 
ual [O II] distributions (on which these predictions rest) from 
our ICBS sample, for each redshift interval. The Dressler- 
Shectman low-redshift cluster and field samples follow the 
"continuous star formation" prediction, but the clusters at 
z 0.4 studied in the Morphs collaboration depart, indicat- 
ing a higher fraction of starbursts (see D04). 
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Fig. 2. Two measures of starburst fraction, SBF, versus redshift for galaxies 
brighter than M*+l. Top: fraction of galaxies with excess EW(H5) > 2.5A; bot- 
tom: fraction of galaxies with SFR > 3 X L,oi X (A//L)/r{universe). Black filled cir- 
cles — ICBS data; open red circles- data from Noeske etal. (2007); red filled cir- 
cles — data from Bell etal. (2005). Black open circles are ICBS data with SFR 
> 10 X L,(ii X {M/L)/T(universe). Using this criterion for a starburst, all three samples 
give consistent results. 

Figure 2 shows that such starbursts make up ^^20-25% of 
starforming galaxies by z ~ 0.6, a much larger fraction than 
today. Moreover, because the "duty cycle" (fraction of the 
time these galaxies are identifiable as starbursts) is almost cer- 
tainly less than 50%, the clear implication is that most star- 
forming galaxies at z> 0.6 have undergone a starburst of at 
least moderate strength. Figure 2 also shows the SBF for a 
factor-of-ten increase rising in parallel. 

Noeske et al. (2007) emphasize the rarity of starbursts over 
a comparable (but somewhat wider) redshift range, basing 
their argument on the narrowness of the SFR/M versus M 
relation. They conclude that no more than one-third of typ- 
ical starforming galaxies have SFR variations greater than a 
factor-of -two, while our sample suggests, after correction for 
the fraction of passive galaxies, that one-third have experi- 
enced a rise of a factor of three or more. Put in these terms, 
the two results appear mildly inconsistent, however, as Figure 
2 illustrates, the Noeske et al. data, when analyzed using the 
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Fig. 3. H(5 versus [O II] equivalent widths for composite spectra, adapted from D04. 
The lines '2dF' and 'CNOC2' show predictions for continuously starforming popula- 
tions, and the points 'Morphs' and 'Dressler-Shectman' cluster and field samples, from 
studies described in D04. 'The low-redshift SWIRE field galaxy sample is the blue point. 
The ICBS field sample measurements are the 5 green points for redshift intervals cen- 
tered at z = 0.25, 0.35, 0.45, 0.55, 0.65. The 5 joined purple points show the prediction 
of EW([0 II] , ) at these 5 redshift intervals for a population of present-epoch star- 
forming galaxies (the SWIRE sample) with very few starbursts, as described in the text. 
The observed ICBS measurements at j = 0.45, 55, &0.65 far exceed these predicted val- 
ues, showing the substantial fraction of starburst galaxies in these field populations. The 
DEEP2 (red) point for field galaxies at ^ ~ 0.83 extends and confirms this result. 

For the new field sample. Figure 3 shows that the z = 
0.25,0.35 slices are little different from the present-epoch 
population, but the z = 0.45,0.55,0.65 slices have much 
stronger H(5 absorption than expected for a population domi- 
nated by continuous star formation. We also can compare our 
EW([0 II] , H(5) values to those from a composite spectrum of 
-1000 field galaxies at z 0.8 from the DEEP2 Galaxy Red- 
shift Survey (Davis et al. 2009), This coadded spectrum was 
produced as described in Weiner et al. (2008) by Jeff New- 
mann and Renbin Yan from galaxies selected with limiting 
Mr comparable to our own sample. The DEEP2 point contin- 
ues the trend of yet greater SBF to even higher redshift, for a 
data set independent of our own, with higher S/N and higher 
spectral resolution. 



In summary, methods based on (1) individual HS strengths, 
(2) the SFR-increase-over-past-average, and (3) (EW([0 II], 
HS ) from composite spectra, all point to an SBF that increases 
markedly with redshift. 

5. DISCUSSION AND CONCLUSIONS 

We have used optical spectra and mid-IR fluxes of field 
galaxies in our ICBS galaxy cluster survey to measure the 
rapid decline of the specific star formation rate, SFR, since 
z ^ 0.7. Although qualitatively consistent with other stud- 
ies, we find that — parameterized as the median specific SFR 

— this decline appears to be steepening toward the present 
epoch. This might be evidence for a kind of "downsizing: " el- 
liptical/bulge formation dominates at 2 < z < 6 and smoothly 
transitions to the building of massive galaxy disks, but as 
this wanes for z < 1, there is insufficient mass in the still- 
starforming dwarf galaxy population to prevent what is es- 
sentially the end of the era of star formation (Dressier 2004). 

We also show that it is not just the SFR, but also its mode 

— the starburst fraction, SBF — that is evolving rapidly since 
z = 0.7. We use three different methods, from discrete and 
composite spectra, to show that the SBF also declines dramat- 
ically over this period. When accounting for the "duty cycle," 
the data suggest that the majority of starforming galaxies at 
z Z 0.5 followed this mode of star formation rather than in the 
steadier, "continuous" mode that dominates today. We em- 
phasize that, unlike short and often-intense nuclear starbursts, 
these are typically of moderate strength — the SFR rising by 
a factor of 3-10. They are likely to be galaxy-wide and of 
longer duration — consistent with the dynamical timescales 
of the larger region involved. 

It seems possible that these two effects — the rapid decline 
in specific SFR and the changing mode of star formation — 
share a common cause. Availability of gas for star forma- 
tion could be the principal agent. If the Schmidt-Kennicutt 
(Kennicutt 1998) law holds, higher SFRs of intermediate- 
redshift galaxies were due to gas fractions higher than the 20- 
40% typical of today's luminous spirals (McCaugh & de Blok 
1997), although how much greater may depend on the rela- 
tive importance of H2 versus HI (Robertson & Kravtsov 2009, 
Krumholz et al. 2009). As discussed by Putnam et al. (2008), 
this greater gas content could be due to initial supply, ac- 
cretion of gas-rich satellites, or resupply from surrounding 
reservoirs of cooling gas, the latter the subject of much re- 
cent discussion (see, e.g., Keres et al. 2009). However, di- 
rect measurements of gas contents for intermediate-redshift 



galaxies are not yet possible, and absorption measurements 
by Prochaska and Wolfe (2009) — albeit of HI alone — sup- 
port a different picture in which gas fractions do not evolve 
strongly and high SFRs are instead supported by high rates of 
gas accretion. 

If higher SFRs were the result of higher gas fractions of 
more than 50% for typical z > 0.5 galaxies, then it is rea- 
sonable to suppose that the higher SBF might also be a re- 
sult. Star formation in gas-rich disks may be unstable if high 
supernovae rates heat disks sufficiently to disrupt conditions 
favorable for star formation. After a gas-cooling time of sev- 
eral hundred million years, rapid star formation could return; 
color-magnitude diagrams for some nearby dwarf galaxies 
show two or more episodes of star formation separated by 
billions of years (e.g, Gallart et al. 1999; Held et al. 2000). 
If resupply by cold gas flows is important to the evolution of 
spiral disks, the higher infall rates of earlier epochs could be 
subject to significant variability, simply a result of the granu- 
larity of the gas density in a complex network of filaments. 

Finally, mergers of gas-rich galaxies might contribute 
to the rising SBF with increasing redshift, although Bell 
et al. (2005) find only a small fraction of their high SFR 
galaxies at z ~ 0.7 in major mergers. Perhaps a higher rate 
of accretion of small satellites at earlier times might have re- 
sulted in more starbursts. Indeed, since these smaller systems 
were probably gas rich, accretion rather than major mergers 
might be the dominant starburst trigger. 

As yet there are no direct means to measure gas fractions 
in z ^ 1 galaxies, but ALMA and the proposed SKA will in 
future provide that capability. Until that time, there is an im- 
portant role for numerical modeling, with its increasing res- 
olution, to explore whether gas fractions, inflow rates, and 
satellite accretion would affect not just the rate of star forma- 
tion but its mode as well. 
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